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Abstract 

The arrival directions of extremely high energy cosmic rays (EHECR) 
above 4 x 10 19 eV, observed by four surface array experiments in the 
northern hemisphere, are examined for coincidences from similar direc- 
tions in the sky. The total number of cosmic rays is 92. A significant 
number of double coincidences (doublet) and triple coincidences (triplet) 
are observed on the supergalactic plane within the experimental angular 
resolution. The chance probability of such multiplets from a uniform dis- 
tribution is less than 1 % if we consider a restricted region within ±10° 
of the supergalactic plane. Though there is still a possibility of chance 
coincidence, the present results on small angle clustering along the su- 
pergalactic plane may be important in interpreting EHECR enigma. An 
independent set of data is required to check our claims. 

PACS Number : 96.40.Pq (Extensive air showers), 98.70.S (Cosmic 
rays(Sources, galactic and extragalactic)), 98.65.-r (Galaxy groups, clus- 
ters, and superclusters; large scale structure of the Universe) 



1 Introduction 



It is well known that extremely high energy cosmic rays (EHECR) are subject to photopion 
production by interactions with primordial photons on traversal through intergalactic space, 
attenuating their energies down to 4xl0 19 eV (40EeV) if their sources are distributed uni- 
formly over cosmological distances. In this case the arrival directions of the highest energy 
cosmic rays around 40EeV are expected to be almost isotropic and some of them above 
40EeV may be correlated to the topological structure of nearby galaxies since the distance 



1 



to their birth places is no further than 50Mpc and the intergalactic magnetic field is at 
most 10~ 9 Gauss If powerful radio galaxies are sources of EHECR's a correlation of 
their arrival directions on the sky with the supergalactic (SG) plane might be expected, 
as suggested by Biermann and Stanev || , since extragalactic radio sources concentrate to- 
wards the SG plane and this concentration extends to at least z ~ 0.02, based on the MRC 
Catalog measured at Molonglo covering the declination band between 18.5° and -18.5° S. 
In 1995 Stanev et al. [Q studied the arrival direction of cosmic rays with energies greater 
than 40EeV and found that they exhibit a good correlation with the direction of the SG 
plane. The magnitude of the observed excess is 2.5^ 2.8(7, in terms of Gaussian probabil- 
ities. They used 42 events of the world data published at that time with 27 events being 
from the Havcrah Park experiment. 

The AG AS A group claimed §| that three doublets of showers with angular separation of less 
than 2.5° (consistent with the experimental resolution) are observed among the 36 events 
above 40EeV, corresponding to a chance probability of 2.9% from a uniform distribution, but 
noted that a significant fraction of EHECR's are uniformly distributed over the observable 
sky. Two of three doublets are observed within 2.0° of the SG plane. Adding the AGASA 
data up to August, 1998 and reevaluating the energies and arrival directions of all AGASA 
events, Takeda et al. (6) found one triplet and three doublets in a total of 47 events. The 
chance probability is smaller than 1 % and the significance increased with the increased 
data set. It is interesting to note that the centroids of the triplet and one of the doublets 
are 0.7 and 0.9 degrees off the SG plane respectively. 

In the northern hemisphere, four surface array experiments at Volcano Ranch (VR), Haverah 
Park (HP), Yakutsk (YK) and Akeno (AG) have been operational so far and, with 81 events 
above 40EeV, we have reported in |tJ that a significant number of doublets and triplet are 
observed around the SG plane with a chance probability less than 1 % supporting the 
AGASA result. 

It is important now to examine whether the significance of such a clustering of EHECR's 
in the sky in correlation with the SG plane increases or not when using the world data set 
which includes the new published AGASA data. If clusters along the SG plane are real, 
source models which do not relate to the SG plane may find difficulty in explaining the 
observations. 

In this report we examine the arrival direction distribution of the events so far published by 
the four experiments in the northern hemisphere. While we report rather low probabilities 
of various occurrences we emphasis that we regard these claims as charting the way for the 
analysis of future, independent data sets rather than as providing conclusive evidence of 
anisotropy. 



2 Analysis and results 



The experimental data reported by VR fill, HP @], YK Jl0| and AG § are used in the 
following analysis. The Fly's Eye events have differing error ellipses, event by event JtJ, 
which makes the estimation of the chance probability complicated and the exposure in 
right ascension of the Fly's Eye instrument is less uniform than that of the ground arrays. 
Hence they are not included in the present analysis. 

Experimental details for the four experiments and conditions for the present analysis are 
summarized in Table 111. Only extensive air showers (EAS) with zenith angles within 45° are 
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used in the present analysis to select EAS of good quality in energy and arrival direction 
determination. The arrival directions are determined by the arrival time of the shower front 
at different detectors, resulting in the estimated directional uncertainties evaluated by each 
experiment as listed in Table [l]for 40EeV events. 

The method of energy estimation is different in each experiment. In this report we use 92 
extensive air showers whose energies are estimated to be above 40EeV by each experiment 
and do not normalize the energy between experiments. 

The observable sky is different for each experiment and is not uniform in declination. In the 
case of the ground array experiments listed in Table |l|, the right ascension distribution is 
almost uniform. The declination distributions for observed events above lOEeV is, however 
different for each experiment as it depends upon the latitude of the array and detector type, 
as shown in Figure 0. Since the triggering efficiency of each experiment becomes uniform 
over their array area above lOEeV, the declination distribution of observed events above 
40EeV in each experiment is similar to that above lOEeV. 



2.1 Galactic and supergalactic latitude distribution 



To study the general arrival direction distribution of 92 events, the latitude distribution 
in galactic (G) and SG coordinates are shown in Figure ||. Solid lines show the sum of 
the expected distributions assuming a uniform incident arrival direction distribution for 
each experiment; a uniform distribution in right ascension and the observed declination 
distributions of Figure [l] are used to obtain these curves. There are no statistically significant 
deviations from uniformity in G coordinates. In the SG coordinate system, there is some 
excess between 0° and 30°, but it is not statistically significant. Most EHECR's, then, are 
found to be uniformly distributed over the observable sky. 



2.2 Clustering of arrival directions of showers 



In Figures |3J and |4], the arrival directions of each event from the four arrays are plotted in 
equatorial and galactic coordinates, respectively. Some coincidences of events coming from 
similar direction are apparent. In the following we estimate the chance probability of such 
coincidences arising from a uniform incident arrival direction distribution. An estimate of 
the space angle scatter due to errors in the arrival direction determination is obtained by 
quadrature addition of the uncertainties of Table ^ and ranges from 2.5° ~ 4.2°, depending 
on the relevant combination of experimental uncertainties. Since 80 % of all events are from 
AGASA and Haverah Park, it is sufficient to examine clusters within 3° ~ 4° for combined 
data set. 

The number of triplets and doublets within space angles of 3°, 4° and 5° are listed in 
Table ||. In counting doublets in this table, each triplet is also decomposed into 3 doublets 
corresponding to the possible pair combinations with differing space angle separations. So 
there are actually 2 triplets and 6 independent doublets within a space angle 3°, but the 
2 triplets are also counted as 6 doublets, making a total of 12 doublets as listed in the 
Table. To estimate the probability of obtaining such clusters by chance from a uniform 
distribution, we simulated the same number of events from each experiment. 1,000,000 data 
sets, each comprising 92 showers, were simulated under the following assumptions : 
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1. The declination distributions of AG, HP and YK events are approximated by smooth 
functions as shown by the solid lines in Figure [j]. The AG declination distribution is 
used to simulate the VR data set. 

2. A uniform distribution in right ascension with the observed declination distributions 
is assumed for all experiments. 

In Figure |5] the frequency distributions of the number of doublets found by simulation of 
1,000,000 data sets within three different space angles separations are shown. The arrow 
shows the observed number of doublets. The chance probability of observing doublets is 
calculated as that fraction of the simulated data sets which have equal, or more doublets 
than observed. 

The chance probabilities of observing triplets or doublets within the three space angles 
separations are listed in Table ||. The probability is small only for clustering within 3° . The 
chance probabilities for clustering within 4° and 5°, are larger than 10 % and hence not 
statistically significant. 



2.3 Correlation with the Supergalactic Plane 



The arrival directions of EHECR of the four ground array experiments are plotted in the 
SG coordinate system in Figure |[ Error circles are also shown; 1.8° in the case of AG, and 
3.0° for HP, VR and YK. 

Note that the two centroid of triplets are in a very narrow region (within 0.9 degrees) 
about the SG plane. An estimate of the chance probabilities of obtaining this number of 
doublets within ±10° from the SG plane is shown in Figure [?]. In this analysis, only the 
number of doublets within ±10° region is counted and compared between observed data 
and simulated data, without taking into account doublets outside this region. In Table |^ 
the chance probabilities of doublets and triplets within ±10° from the SG plane are listed. 
Whilst the choice of a cut at 10 degrees about the SG plane is arbitrary, we note that the 
chance probabilities for doublets within all space angle separations, and for triplets within 
the two smaller space angle cuts, are intriguingly below 1 %. 



3 Discussion 



As Figure [2] indicates, there is about a 20 % excess of events in the latitude band ~ 30° 
about the SG plane. Though this excess is only a one a effect, it may contribute to the 
apparently significant clustering estimated above, so we have also calculated the chance 
probability of multiplets, given a uniform distribution within this band but with a total event 
number excess of 20 %. The result is listed in Table ||. As expected the statistical evidence 
for clustering is reduced somewhat (the chance probability is raised by approximately a 
factor 2), but for doublets within 4° the chance probability remains below 1 %. 

To examine the effect of the choice of coordinates, the analysis is repeated for the G co- 
ordinate system. There are no multiplets found, even for the largest space angles (within 
5°) and within ±10° of the G plane. The average number of doublets expected by chance 
is 2.8. The Poisson probability of observing no multiplets is thus 6 %. Considering that the 
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number of events within ±10° of the G plane is about 30 % less than that expected from 
uniformity (see Fig. 2), this probability is statistically reasonable. 



The details of each event which are members of one of the clusters within 4° are listed in 
Table ^. Here one HP event is counted twice (doublet #5 and #8), since it forms a doublet 
with another HP event, and with a VR event. Clearly some, at least, of these multiplets must 
arise by chance. In an attempt to select genuinely coincident events from sources at limited 
distances, it may be better to only use events above a slightly higher energy threshold, to 
avoid the systematic differences in energy determination of each experiment. Selecting only 
events above 5 x 10 19 eV, two triplets and one doublet (#7) remain, from a total number of 
51 events. The chance probability of the triplet is now 1.1 %, falling to 0.1 % if we examine 
only the restricted region within ±10° of the SG plane. 



The order of arrival times with respect to energy of the events in triplct^l are not as 
expected from a recent, nearby bursting source where cosmic rays are accerelated in a 
short time According to simulations on the propagation of protons through both the 
intergalactic and Galactic magnetic fields by Medina Tanco jl2| , the extragalactic magnetic 
field must be much smaller than the present upper limit of 10~ 9 G to explain the present 
clustering within experimental angular resolution. 

In the case of triplet #2, the direction is consistent with the Ursa-Major II cluster of 
galaxies. The magnetic field strength in this cluster of galaxies is possibly of the order of 
sub-/zG, similar to that observed in the Coma cluster If so, each member of the triplet 
might be a gamma-ray, because protons may not be collimated within 2.5° and the mean 
distance of travel before decay for neutrons is 1 Mpc for 10 20 eV. Since the cross-section 
of pair creation by gamma-rays of this energy and the probability of bremsstrahlung of the 
resulting electrons are suppressed due to the relativistic contraction of the atmosphere (LPM 
effect [|l5| [l6|]), the longitudinal development of a gamma-ray shower is greatly depressed 
and delayed. A ground array, therefore, might grossly underestimate the primary energy of 
such a shower if it is estimated by the local particle density around 600m from the shower 
core |(J. However, if the geomagnetic field component normal to the arrival direction is 
large enough, electron-positron pair creation in the geomagnetic field occurs far from the 
earth and a cascade develops in the stratoshpereJlTj. Therefore in the northern hemisphere 
we might observe gamma-ray showers in the highest energy region mainly from a northerly 
direction jL8|. The terrestrial arrival directions of the three members of triplet jj=2 are 
indeed from the north, which is at least consistent with this conjecture. However, it should 
be noted that the attenuation length of gamma-rays of energies above 4 x 10 19 eV due to 
the interaction on radio photons is less than 10 Mpc |l3 and hence a source distance limit 
also applies for gamma-rays. 

If each member of the triplets are protons, coming from the same source, then the inter- 
vening magnetic field must be so weak that the particles are hardly deflected, or there is 
magnetic focusing in the magnetic field structure of the Local Supercluster as demonstrated 
by Lemoine et al. Jl9|] . 

In Figure ||, the directions of triplets (open squares) and doublets (open circles) are plotted 
on the CfA galaxy ^TJ| distribution within 100 Mpc. There are no multiplets from the most 
crowded region (Virgo Cluster) and there seems no correlation with the density of nearby 
galaxies. In the following we look for any source candidate for the triplets #1 and #2, and 
the doublet #7. 
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The AG highest energy event and HP 10 20 eV event are members of the triplet #1 and this 
triplet may be related to a nearby source. Mrk 359 (1= 134.60°, b=-42.87°) with z=0.017 
(68 Mpc assuming Ho=75 km/s/Mpc) is within 2.3° from the direction of the centroid of 
this triplet. This direction is also within 3 and 6 degrees of a clusterings of events with 
energies above 10 19 eV claimed by Chi et al. (1~ 133°, b~ -40°) [gl] and Efimov and 
Mikhailov (RA=27°, Dec=18° or 1=143°, b=-42°) ||. Al-Dargazelli et al. [|| pointed 
that the colliding galaxy pair VV338 (1=138°, b=-34°) (N672 and U1249) may be related 
to the clustering around the region 1=135° and b=-35°, where they identified an apparent 
clustering above 10 19 eV. Though VV338 is very close, 5.7 Mpc, the angular separation 
from triplet #1 is large (9.0°). 

Triplet #2 comprises three AG events of similar energies, (5-8) xlO 19 eV, and is within 
2.4° from the direction of the interacting galaxy, VV144 (Mkn40 or Arpl51; 1=147.03°, 
b=58.54°). This is a Seyfert galaxy with z=0.020, corresponding to about 81 Mpc. Triplet 
#2 is on top of a maximum in the arrival direction probability simulated by Medina Tanco 
[fj"2| for sources located between 20 and 50 Mpc. These simulations assume that the luminous 
matter in the nearby universe tracks the distribution of cosmic ray sources and modulates 
the intensity of IGMF. 

In the direction of doublet #7, there are interacting galaxies VV101 (1=87.06°, b=33.75°, 
100 Mpc) with a space angle separation of 1.3° and VV89 (1=88.14°, b=35.51°, 14.5 Mpc) 
with a separation of 2.7°. 

From the above discussion, if the extragalactic magnetic field is much weaker than the 
present upper limit of 10~ 9 G, then triplets #1 and #2 and doublet #7 are possibly corre- 
lated with interacting galaxies, as pointed out by Al-Dargazelli et al ^3|. Clearly, the next 
generation experiments are needed to confirm such an association. 

Recently Lemoine et al. made extensive simulations of the propagation of protons in 
a large scale magnetic field of strength ~ 0.05 — 0.5/xG in the Local Supercluster. They 
found an 8-20 % probability of detecting 5 doublets above 40 EeV in the AGASA data set 
of 47 events, due to magnetic focusing. It is important to examine what kind of conditions 
are required for the magnetic field configuration in the Local Supercluster to explain the 
present results of two triplets along the SG plane. 



4 Conclusion 



A number of collimated triplets and doublets within 4° (the approximate angular resolution 
of combined data set) are observed. The chance probability is of the order of 10 % and is 
not statistically significant. However, the chance of observing triplets and doublets within 
±10° of the SG plane is less than 1 %. It should be noted that the probability of observing 
the multiplets above 5xl0 19 eV is about 1 % and it is less than 0.1 % if attention is limited 
to within ±10° of the SG plane. Though there is still a possibility of chance coincidence, 
we should pay more attention on the small angle clustering along the SG plane in the 
interpretation of EHECR enigma. 

We expect detailed investigation with much better angular resolution and much higher 
statistics from future experiments. In particular we note that the probabilities given in this 
paper are 'a postiori' in that hints of clustering have been reported in earlier analysis and 
these data have been included here. There are insufficient events for an independent test 
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of our observations. An independent set awaits new projects such as Hi-Res, Auger and 
Telescope Array. 
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Fig.l Declination distribution of showers of energies above lOEcV and zenith angles smaller 

than 45° of AG, YK and HP. 
Fig. 2 Galactic (left) and supergalactic latitude (right) distribution of arrival directions of 

92 cosmic rays from four ground array experiments. Energy thresholds are 40EcV. 

Solid lines are the sum of expected distributions of each experiment for a uniform 

distribution. 

Fig. 3 Sky map of 92 events above 40EeV in equatorial coordinates. Squares - AGASA, 
Triangles - Haverah Park, Circles - Yakutsk and Stars - Volcano Ranch. The region 
of sky observable in each experiment is shown by dashed lines from the top YK, HP, 
AG and VR. 

Fig. 4 Sky map of 92 events above 40EeV in galactic coordinates. The symbols are the same 
as used in figure 3. 

Fig. 5 Frequency distribution of the number of doublets in 1,000,000 simulated data sets 
within several space angles. The arrow mark shows the observed number of doublets. 
Chance probability of the observed number of doublets is calculated by summing the 
number of data sets in the hatched region. 

Fig. 6 Arrival directions with error circles of EHECR in supergalactic coordinates. The 
total number of events is 92. 

Fig. 7 Probability of observing doublets and triplets within 10° from the SG plane with 
1,000,000 simulated data sets for 92 events. 

Fig. 8 The location of galaxies within 100 Mpc in G coordinates from the CfA 1995 cat- 
alogue and the clusters within 4°. The squares show triplets and the circles show 
doublets. 
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Experiment 


Longitude 


Latitude 


Number of 
events 


Zenith 
angle 


Error in arrival 
determination 


AGASA 


138°30'E 


38°47'N 


47 


< 45° 


1.8° 


Havcrah Park 


1°38'W 


53°58'N 


27 


< 45° 


3° 


Yakutsk 


129°24'E 


61°42'N 


12 


< 45° 


3° 


Volcano Ranch 


106°47'W 


35°09'N 


6 


< 45° 


(3°) 


Total 






92 







Table 1: Experimental sites and number of events above J^OEeV. 



Combination 


Resolution (degrees) 


AG-AG 


2.5 


AG-YK, AG-HP, AG-VR 


3.5 


HP-HP, YK-YK, VR-VR 


4.2 


HP-YK, HP-VR, YK-VR 


4.2 



Table 2: Combined space angle resolution from angular resolution of each experiment. 





Doublet 


Triplet 


space angle 


observed 


probability 


observed 


probability 




number 




number 




< 3.0° 


12 


1.5% 


2 


1.4% 


< 4.0° 


14 


13.4% 


2 


8.3% 


< 5.0° 


20 


15.9% 


3 


11.8% 



Table 3: Probability of observing doublets and triplets with 1,000,000 simulated data sets. 
Each triplet is counted as three doublets. 





Doublet 


Triplet 


space angle 


observed 
number 


probability 


observed 
number 


probability 


< 3.0° 


8 


0.1% 


2 


0.2% 


< 4.0° 


9 


0.3% 


2 


0.9% 


< 5.0° 


11 


0.6% 


2 


3.1% 



Table 4: Probability of observing doublets and triplets within ±10° from the SG plane for 
92 events. 





Doublet 


Triplet 


space angle 


observed 
number 


probability 


observed 
number 


probability 


< 3.0° 


8 


0.2% 


2 


0.3% 


< 4.0° 


9 


0.8% 


2 


1.7% 


< 5.0° 


11 


1.5% 


2 


5.4% 



Table 5: Probability of observing doublets and triplets when an excess of 20 % in the 0° ~ 30° 
range of SG latitide is assumed. 
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V_/l LLo Lt-I 






Log E 




JJct. 


l 


u 


O. VJ.ljIlg. 


O. Vjr.ljcLL. 


Triplet #1 


HP 


810105 


19.99 


20.00 


20.00 


132.70 


-41.70 


318.10 


-0.79 




AG 


C\ O 1 O C\ o 

931203 


20.33 


18.91 


21.07 


130.48 


-41.44 


318.11 


0.89 




AG 


951029 


19.71 


18.53 


20.03 


130.18 


-42.51 


317.02 


0.93 


Triplet #2 


AG 


AOAOA1 

920801 


19.74 


172.30 


57.14 


143.20 


56.65 


56.82 


2.04 




AG 


950126 


19.89 


168.65 


57.58 


145.53 


55.10 


55.51 


0.51 




A f ' 

AG 


C\OC\ Af\ A 

980404 


i n to 

19.73 


1 CO A A 

168.44 


55.99 


147.51 


56.23 


£ G OA 

56.84 


A O T 

-0.37 


Doublet #1 


AG 


910420 


19.64 


284.90 


47.79 


77.88 


18.45 


24.95 


57.83 




AG 


940706 


20.03 


281.36 


48.32 


77.58 


20.86 


29.35 


57.26 


Doublet #2 


AG 


860105 


19.74 


69.03 


30.15 


170.08 


-11.50 


350.38 


-33.33 




AG 


951115 


19.69 


70.39 


29.85 


171.09 


-10.79 


351.23 


-34.31 


Doublet #3 


HP 


860315 


19.71 


267.00 


77.00 


108.50 


30.10 


30.83 


27.99 




AG 


960513 


19.68 


269.05 


74.12 


105.11 


29.79 


31.09 


30.94 


Doublet #4 


HP 


720525 


19.65 


239.00 


79.00 


113.30 


34.60 


35.05 


23.27 




YK 


911201 


19.62 


235.40 


79.80 


114.60 


34.60 


34.88 


22.22 


Doublet #5 


VR 


610319 


19.73 


154.10 


66.70 


143.00 


44.30 


44.59 


0.35 




HP 


850313 


19.62 


157.00 


65.00 


143.60 


46.30 


46.63 


0.24 


Doublet #6 


HP 


661008 


19.67 


164.00 


50.00 


159.00 


58.80 


61.08 


-5.53 




YK 


750317 


19.67 


163.70 


52.90 


154.90 


56.80 


58.45 


-4.16 


Doublet #7 


HP 


740228 


19.86 


264.00 


58.00 


86.36 


32.52 


41.02 


45.22 




AG 


980330 


19.84 


259.16 


56.32 


84.39 


35.17 


45.44 


45.35 


Doublet #8 


HP 


760206 


19.62 


165.00 


64.00 


140.98 


49.43 


49.49 


2.41 




HP 


850313 


19.62 


157.00 


65.00 


143.60 


46.30 


46.63 


0.24 



Table 6: Event lists of members of clusters within 4° from 4 surface experiments 
AGASA(AG), Haverah Park(HP), Volcano Ranch(VR) and Yaktuk(YK). 



Figure 1: Declination distribution of showers of energies above lOEeV and zenith angles 
smaller than 45° of AG, YK and HP. 
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Figure 2: Galactic (left) and supergalactic latitude (right) distribution of arrival directions 
of 92 cosmic rays from four ground array experiments. Energy thresholds are J^OEeV. Solid 
lines are the sum of expected distributions of each experiment for a uniform distribution. 



Figure 3: Sky map of 92 events above J^OEeV in equatorial coordinates. Squares - AGASA, 
Triangles - Haverah Park, Circles - Yakutsk and Stars - Volcano Ranch. The region of sky 
observable in each experiment is shown by dashed lines from the top YK, HP, AG and VR. 
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Figure 4: Sky map of 92 events above J^OEeV in galactic coordinates. The symbols are the 
same as used in figure 3. 



Figure 5: Frequency distribution of the number of doublets in 1,000,000 simulated data sets 
within several space angles. The arrow mark shows the observed number of doublets. Chance 
probability of the observed number of doublets is calculated by summing the number of data 
sets in the hatched region. 
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Figure 6: Arrival directions with error circles of EHECR in supergalactic coordinates. The 
total number of events is 92. 



Figure 7: Probability of observing doublets and triplets within 10° from the SG plane with 
1,000,000 simulated data sets for 92 events. 
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Figure 8: The location of galaxies within 100 Mpc in G coordinates from the CfA 1995 
catalogue and the clusters within 4°. The squares show triplets and the circles show doublets. 
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Abstract 

The arrival directions of extremely high energy cosmic rays (EHECR) above 
4 X 10 19 eV, observed by four surface array experiments in the northern hemisphere, 
are examined for coincidences from similar directions in the sky. The total number 
of cosmic rays is 92. A significant number of double coincidences (doublet) and triple 
coincidences (triplet) are observed on the supergalactic plane within the experimental 
angular resolution. The chance probability of such multiplets from a uniform distribu- 
tion is less than 1 % if we consider a restricted region within ±10° of the supergalactic 
plane. Though there is still a possibility of chance coincidence, the present results on 
small angle clustering along the supergalactic plane may be important in interpreting 
EHECR enigma. An independent set of data is required to check our claims. 

PACS Number : 96.40.Pq (Extensive air showers), 98.70.S (Cosmic rays(Sources, 
galactic and extragalactic)), 98.65.-r (Galaxy groups, clusters, and superclusters; large 
scale structure of the Universe) 

1 Introduction 

It is well known that extremely high energy cosmic rays (EHECR) are subject to photopi- 
on production by interactions with primordial photons on traversal through intergalactic 
space, attenuating their energies down to 4xl0 19 eV (40EeV) if their sources are distribut- 
ed uniformly over cosmological distances. In this case the arrival directions of the highest 
energy cosmic rays around 40EeV are expected to be almost isotropic and some of them 
above 40EeV may be correlated to the topological structure of nearby galaxies since the 
distance to their birth places is no further than 50Mpc and the intergalactic magnetic 
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field is at most 10~ 9 Gauss [1]. If powerful radio galaxies are sources of EHECR's a cor- 
relation of their arrival directions on the sky with the supergalactic (SG) plane might 
be expected, as suggested by Biermann and Stanev [2], since extragalactic radio sources 
concentrate towards the SG plane and this concentration extends to at least z ~ 0.02, 
based on the MRC Catalog measured at Molonglo covering the declination band between 
18.5° and -18.5° [3]. In 1995 Stanev et al. [4] studied the arrival direction of cosmic rays 
with energies greater than 40EeV and found that they exhibit a good correlation with the 
direction of the SG plane. The magnitude of the observed excess is 2.5~ 2.8<r, in terms 
of Gaussian probabilities. They used 42 events of the world data published at that time 
with 27 events being from the Haverah Park experiment. 

The AGASA group claimed [5] that three doublets of showers with angular separation 
of less than 2.5° (consistent with the experimental resolution) are observed among the 
36 events above 40EeV, corresponding to a chance probability of 2.9% from a uniform 
distribution, but noted that a significant fraction of EHECR's are uniformly distributed 
over the observable sky. Two of three doublets are observed within 2.0° of the SG plane. 
Adding the AGASA data up to August, 1998 and reevaluating the energies and arrival 
directions of all AGASA events, Takeda et al. [6] found one triplet and three doublets 
in a total of 47 events. The chance probability is smaller than 1 % and the significance 
increased with the increased data set. It is interesting to note that the centroids of the 
triplet and one of the doublets are 0.7 and 0.9 degrees off the SG plane respectively. 

In the northern hemisphere, four surface array experiments at Volcano Ranch (VR), Hav- 
erah Park (HP), Yakutsk (YK) and Akeno (AG) have been operational so far and, with 
81 events above 40EeV, we have reported in [7] that a significant number of doublets 
and triplet are observed around the SG plane with a chance probability less than 1 % 
supporting the AGASA result. 

It is important now to examine whether the significance of such a clustering of EHECR's 
in the sky in correlation with the SG plane increases or not when using the world data set 
which includes the new published AGASA data. If clusters along the SG plane are real, 
source models which do not relate to the SG plane may find difficulty in explaining the 
observations. 

In this report we examine the arrival direction distribution of the events so far published by 
the four experiments in the northern hemisphere. While we report rather low probabilities 
of various occurrences we emphasis that we regard these claims as charting the way for 
the analysis of future, independent data sets rather than as providing conclusive evidence 
of anisotropy. 
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2 Analysis and results 



The experimental data reported by VR [11], HP [4], YK [10] and AG [6] are used in the 
following analysis. The Fly's Eye events have differing error ellipses, event by event [7], 
which makes the estimation of the chance probability complicated and the exposure in 
right ascension of the Fly's Eye instrument is less uniform than that of the ground arrays. 
Hence they are not included in the present analysis. 

Experimental details for the four experiments and conditions for the present analysis are 
summarized in Table 1. Only extensive air showers (EAS) with zenith angles within 45° are 
used in the present analysis to select EAS of good quality in energy and arrival direction 
determination. The arrival directions are determined by the arrival time of the shower 
front at different detectors, resulting in the estimated directional uncertainties evaluated 
by each experiment as listed in Table 1 for 40EeV events. 

The method of energy estimation is different in each experiment. In this report we use 92 
extensive air showers whose energies are estimated to be above 40EeV by each experiment 
and do not normalize the energy between experiments. 

The observable sky is different for each experiment and is not uniform in declination. In the 
case of the ground array experiments listed in Table 1, the right ascension distribution is 
almost uniform. The declination distributions for observed events above lOEeV is, however 
different for each experiment as it depends upon the latitude of the array and detector 
type, as shown in Figure 1. Since the triggering efficiency of each experiment becomes 
uniform over their array area above lOEeV, the declination distribution of observed events 
above 40EeV in each experiment is similar to that above lOEeV. 

2.1 Galactic and supergalactic latitude distribution 

To study the general arrival direction distribution of 92 events, the latitude distribution 
in galactic (G) and SG coordinates are shown in Figure 2. Solid lines show the sum of 
the expected distributions assuming a uniform incident arrival direction distribution for 
each experiment; a uniform distribution in right ascension and the observed declination 
distributions of Figure 1 are used to obtain these curves. There are no statistically sig- 
nificant deviations from uniformity in G coordinates. In the SG coordinate system, there 
is some excess between 0° and 30°, but it is not statistically significant. Most EHECR's, 
then, are found to be uniformly distributed over the observable sky. 
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2.2 Clustering of arrival directions of showers 

In Figures 3 and 4, the arrival directions of each event from the four arrays are plotted in 
equatorial and galactic coordinates, respectively. Some coincidences of events coming from 
similar direction are apparent. In the following we estimate the chance probability of such 
coincidences arising from a uniform incident arrival direction distribution. An estimate of 
the space angle scatter due to errors in the arrival direction determination is obtained by 
quadrature addition of the uncertainties of Table 2 and ranges from 2.5° ~ 4.2°, depending 
on the relevant combination of experimental uncertainties. Since 80 % of all events are 
from AGASA and Haverah Park, it is sufficient to examine clusters within 3° ~ 4° for 
combined data set. 

The number of triplets and doublets within space angles of 3°, 4° and 5° are listed in 
Table 3. In counting doublets in this table, each triplet is also decomposed into 3 doublets 
corresponding to the possible pair combinations with differing space angle separations. So 
there are actually 2 triplets and 6 independent doublets within a space angle 3°, but the 
2 triplets are also counted as 6 doublets, making a total of 12 doublets as listed in the 
Table. To estimate the probability of obtaining such clusters by chance from a uniform 
distribution, we simulated the same number of events from each experiment. 1,000,000 
data sets, each comprising 92 showers, were simulated under the following assumptions : 

1. The declination distributions of AG, HP and YK events are approximated by smooth 
functions as shown by the solid lines in Figure 1. The AG declination distribution is 
used to simulate the VR data set. 

2. A uniform distribution in right ascension with the observed declination distributions 
is assumed for all experiments. 

In Figure 5 the frequency distributions of the number of doublets found by simulation of 
1,000,000 data sets within three different space angles separations are shown. The arrow 
shows the observed number of doublets. The chance probability of observing doublets is 
calculated as that fraction of the simulated data sets which have equal, or more doublets 
than observed. 

The chance probabilities of observing triplets or doublets within the three space angles 
separations are listed in Table 3. The probability is small only for clustering within 3°. 
The chance probabilities for clustering within 4° and 5°, are larger than 10 % and hence 
not statistically significant. 
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2.3 Correlation with the Supergalactic Plane 

The arrival directions of EHECR of the four ground array experiments are plotted in the 
SG coordinate system in Figure 6. Error circles are also shown; 1.8° in the case of AG, 
and 3.0° for HP, VR and YK. 

Note that the two centroid of triplets are in a very narrow region (within 0.9 degrees) 
about the SG plane. An estimate of the chance probabilities of obtaining this number of 
doublets within ±10° from the SG plane is shown in Figure 7. In this analysis, only the 
number of doublets within ±10° region is counted and compared between observed data 
and simulated data, without taking into account doublets outside this region. In Table 4 
the chance probabilities of doublets and triplets within ±10° from the SG plane are listed. 
Whilst the choice of a cut at 10 degrees about the SG plane is arbitrary, we note that 
the chance probabilities for doublets within all space angle separations, and for triplets 
within the two smaller space angle cuts, are intriguingly below 1 %. 

3 Discussion 

As Figure 2 indicates, there is about a 20 % excess of events in the latitude band ~ 
30° about the SG plane. Though this excess is only a one a effect, it may contribute 
to the apparently significant clustering estimated above, so we have also calculated the 
chance probability of multiplets, given a uniform distribution within this band but with 
a total event number excess of 20 %. The result is listed in Table 5. As expected the 
statistical evidence for clustering is reduced somewhat (the chance probability is raised 
by approximately a factor 2), but for doublets within 4° the chance probability remains 
below 1 %. 

To examine the effect of the choice of coordinates, the analysis is repeated for the G 
coordinate system. There are no multiplets found, even for the largest space angles (within 
5°) and within ±10° of the G plane. The average number of doublets expected by chance 
is 2.8. The Poisson probability of observing no multiplets is thus 6 %. Considering that 
the number of events within ±10° of the G plane is about 30 % less than that expected 
from uniformity (see Fig. 2), this probability is statistically reasonable. 

The details of each event which are members of one of the clusters within 4° are listed in 
Table 6. Here one HP event is counted twice (doublet ^5 and #8), since it forms a doublet 
with another HP event, and with a VR event. Clearly some, at least, of these multiplets 
must arise by chance. In an attempt to select genuinely coincident events from sources 
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at limited distances, it may be better to only use events above a slightly higher energy 
threshold, to avoid the systematic differences in energy determination of each experiment. 
Selecting only events above 5 X 10 19 eV, two triplets and one doublet (#7) remain, from 
a total number of 51 events. The chance probability of the triplet is now 1.1 %, falling to 
0.1 % if we examine only the restricted region within ±10° of the SG plane. 

The order of arrival times with respect to energy of the events in triplet^l are not as 
expected from a recent, nearby bursting source where cosmic rays are accerelated in a 
short time [8]. According to simulations on the propagation of protons through both 
the intergalactic and Galactic magnetic fields by Medina Tanco [12], the extragalactic 
magnetic field must be much smaller than the present upper limit of 10~ 9 G to explain 
the present clustering within experimental angular resolution. 

In the case of triplet $2, the direction is consistent with the Ursa-Major II cluster of 
galaxies. The magnetic field strength in this cluster of galaxies is possibly of the order of 
sub-//G, similar to that observed in the Coma cluster [13]. If so, each member of the triplet 
might be a gamma-ray, because protons may not be collimated within 2.5° and the mean 
distance of travel before decay for neutrons is 1 Mpc for 10 20 eV. Since the cross-section of 
pair creation by gamma-rays of this energy and the probability of bremsstrahlung of the 
resulting electrons are suppressed due to the relativistic contraction of the atmosphere 
(LPM effect [15, 16]), the longitudinal development of a gamma-ray shower is greatly 
depressed and delayed. A ground array, therefore, might grossly underestimate the primary 
energy of such a shower if it is estimated by the local particle density around 600m from 
the shower core [9, 6]. However, if the geomagnetic field component normal to the arrival 
direction is large enough, electron-positron pair creation in the geomagnetic field occurs far 
from the earth and a cascade develops in the stratoshpere[17]. Therefore in the northern 
hemisphere we might observe gamma-ray showers in the highest energy region mainly 
from a northerly direction [18]. The terrestrial arrival directions of the three members of 
triplet ^2 are indeed from the north, which is at least consistent with this conjecture. 
However, it should be noted that the attenuation length of gamma-rays of energies above 
4 X 10 19 eV due to the interaction on radio photons is less than 10 Mpc [14] and hence a 
source distance limit also applies for gamma-rays. 

If each member of the triplets are protons, coming from the same source, then the inter- 
vening magnetic field must be so weak that the particles are hardly deflected, or there 
is magnetic focusing in the magnetic field structure of the Local Supercluster as demon- 
strated by Lemoine et al. [19]. 

In Figure 8, the directions of triplets (open squares) and doublets (open circles) are plotted 
on the CfA galaxy [20] distribution within 100 Mpc. There are no multiplets from the 
most crowded region (Virgo Cluster) and there seems no correlation with the density of 
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nearby galaxies. In the following we look for any source candidate for the triplets $1 and 
#2, and the doublet #7. 

The AG highest energy event and HP 10 20 eV event are members of the triplet ftl and this 
triplet may be related to a nearby source. Mrk 359 (1= 134.60°, b=-42.87°) with z=0.017 
(68 Mpc assuming Hq=75 km/s/Mpc) is within 2.3° from the direction of the centroid of 
this triplet. This direction is also within 3 and 6 degrees of a clusterings of events with 
energies above 10 19 eV claimed by Chi et al. (1~ 133°, b~ -40°) [21] and Efimov and 
Mikhailov (RA=27°, Dec=18° or 1=143°, b=-42°) [22]. Al-Dargazelli et al. [23] pointed 
that the colliding galaxy pair VV338 (1=138°, b=-34°) (N672 and U1249) may be related 
to the clustering around the region 1=135° and b=-35°, where they identified an apparent 
clustering above 10 19 eV. Though VV338 is very close, 5.7 Mpc, the angular separation 
from triplet #1 is large (9.0°). 

Triplet ^2 comprises three AG events of similar energies, (5-8) xlO 19 eV, and is within 
2.4° from the direction of the interacting galaxy, VV144 (Mkn40 or Arpl51; 1=147.03°, 
b=58.54°). This is a Seyfert galaxy with z=0.020, corresponding to about 81 Mpc. Triplet 
^2 is on top of a maximum in the arrival direction probability simulated by Medina 
Tanco [12] for sources located between 20 and 50 Mpc. These simulations assume that 
the luminous matter in the nearby universe tracks the distribution of cosmic ray sources 
and modulates the intensity of IGMF. 

In the direction of doublet $7, there are interacting galaxies VV101 (1=87.06°, b=33.75°, 
100 Mpc) with a space angle separation of 1.3° and VV89 (1=88.14°, b=35.51°, 14.5 Mpc) 
with a separation of 2.7°. 

From the above discussion, if the extragalactic magnetic field is much weaker than the 
present upper limit of 10~ 9 G, then triplets ftl and ^2 and doublet ^7 are possibly 
correlated with interacting galaxies, as pointed out by Al-Dargazelli et al [23]. Clearly, 
the next generation experiments are needed to confirm such an association. 

Recently Lemoine et al. [19] made extensive simulations of the propagation of protons in 
a large scale magnetic field of strength ~ 0.05 — 0.5//G in the Local Supercluster. They 
found an 8-20 % probability of detecting 5 doublets above 40 EeV in the AGASA data set 
of 47 events, due to magnetic focusing. It is important to examine what kind of conditions 
are required for the magnetic field configuration in the Local Supercluster to explain the 
present results of two triplets along the SG plane. 
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4 Conclusion 



A number of collimated triplets and doublets within 4° (the approximate angular resolu- 
tion of combined data set) are observed. The chance probability is of the order of 10 % 
and is not statistically significant. However, the chance of observing triplets and doublets 
within ±10° of the SG plane is less than 1 %. It should be noted that the probability 
of observing the multiplets above 5xl0 19 eV is about 1 % and it is less than 0.1 % if 
attention is limited to within ±10° of the SG plane. Though there is still a possibility of 
chance coincidence, we should pay more attention on the small angle clustering along the 
SG plane in the interpretation of EHECR enigma. 

We expect detailed investigation with much better angular resolution and much higher 
statistics from future experiments. In particular we note that the probabilities given in 
this paper are 'a postiori' in that hints of clustering have been reported in earlier analysis 
and these data have been included here. There are insufficient events for an independent 
test of our observations. An independent set awaits new projects such as Hi-Res, Auger 
and Telescope Array. 
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Fig.l Declination distribution of showers of energies above lOEeV and zenith angles 

smaller than 45° of AG, YK and HP. 
Fig. 2 Galactic (left) and supergalactic latitude (right) distribution of arrival directions 

of 92 cosmic rays from four ground array experiments. Energy thresholds are 40EeV. 

Solid lines are the sum of expected distributions of each experiment for a uniform 

distribution. 

Fig. 3 Sky map of 92 events above 40EeV in equatorial coordinates. Squares - AGASA, 
Triangles - Haverah Park, Circles - Yakutsk and Stars - Volcano Ranch. The region 
of sky observable in each experiment is shown by dashed lines from the top YK, HP, 
AG and VR. 

Fig. 4 Sky map of 92 events above 40EeV in galactic coordinates. The symbols are the 

same as used in figure 3. 
Fig. 5 Frequency distribution of the number of doublets in 1,000,000 simulated data sets 

within several space angles. The arrow mark shows the observed number of doublets. 

Chance probability of the observed number of doublets is calculated by summing the 

number of data sets in the hatched region. 
Fig. 6 Arrival directions with error circles of EHECR in supergalactic coordinates. The 

total number of events is 92. 
Fig. 7 Probability of observing doublets and triplets within 10° from the SG plane with 

1,000,000 simulated data sets for 92 events. 
Fig. 8 The location of galaxies within 100 Mpc in G coordinates from the CfA 1995 

catalogue and the clusters within 4°. The squares show triplets and the circles show 

doublets. 
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Experiment 


Longitude 


Latitude 


Number of 
events 


Zenith 
angle 

O 


Error in arrival 
determination 


AGASA 


138°30'E 


38°47'N 


47 


< 45° 


1.8° 


Haver ah Park 


1°38'W 


53°58'N 


27 


< 45° 


3° 


Yakutsk 


129°24'E 


61°42'N 


12 


< 45° 


3° 


Volcano Ranch 


106°47'W 


35°09'N 


6 


< 45° 


(3°) 


Total 






92 







Table 1: Experimental sites and number of events above J^OEeV. 



Combination 


Resolution (degrees) 


AG-AG 


2.5 


AG-YK, AG-HP, AG-VR 


3.5 


HP-HP, YK-YK, VR-VR 


4.2 


HP-YK, HP-VR, YK-VR 


4.2 



Table 2: Combined space angle resolution from angular resolution of each experiment. 





Doublet 


Triplet 


space angle 


observed 
number 


probability 


observed 
number 


probability 


< 3.0° 


12 


1.5% 


2 


1.4% 


< 4.0° 


14 


13.4% 


2 


8.3% 


< 5.0° 


20 


15.9% 


3 


11.8% 



Table 3: Probability of observing doublets and triplets with 1,000,000 simulated data sets. 
Each triplet is counted as three doublets. 





Doublet 


Triplet 


space angle 


observed 
number 


probability 


observed 
number 


probability 


< 3.0° 


8 


0.1% 


2 


0.2% 


< 4.0° 


9 


0.3% 


2 


0.9% 


< 5.0° 


11 


0.6% 


2 


3.1% 



Table 4: Probability of observing doublets and triplets within ±10° from the SG plane for 
92 events. 





Doublet 


Triplet 


space angle 


observed 
number 


probability 


observed 
number 


probability 


< 3.0° 


8 


0.2% 


2 


0.3% 


< 4.0° 


9 


0.8% 


2 


1.7% 


< 5.0° 


11 


1.5% 


2 


5.4% 



Table 5: Probability of observing doublets and triplets when an excess of 20 % in the 
0° ~ 30° range of SG latitide is assumed. 
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f lief 

VjlUo lei 




JJalc 


T r>a F 
IjOti Hi 


R A 




1 


u 

u 


o.vj.ljllg. 


S G T 


Triplet #1 


HP 


810105 


19.99 


20.00 


20.00 


132.70 


-41.70 


318.10 


-0.79 




AG 


931203 


20.33 


18.91 


21.07 


130.48 


At A A 

-41.44 


318.11 


0.89 




AG 


951029 


19.71 


18.53 


20.03 


130.18 


-42.51 


317.02 


0.93 


Iriplet #2 


AG 


920801 


19.74 


172.30 


57.14 


143.20 


56.65 


56.82 


d n a 

2.04 




AG 


950126 


19.89 


168.65 


57.58 


145.53 


55.10 


55.51 


0.51 




AG 


f\0 C\ A C\ A 

980404 


1 fl TO 

19.73 


1 C O A A 

lbo.44 


55.99 


1 A T C 1 

147.51 


56.23 


C £ OA 

56.84 


-0.37 


Doublet #1 


AG 


910420 


19.64 


284.90 


47.79 


77.88 


18.45 


24.95 


57.83 




AG 


940706 


20.03 


281.36 


48.32 


77.58 


20.86 


29.35 


57.26 


Doublet #2 


AG 


860105 


19.74 


69.03 


30.15 


170.08 


-11.50 


350.38 


-33.33 




AG 


951115 


19.69 


70.39 


29.85 


171.09 


-10.79 


351.23 


-34.31 


Doublet #3 


HP 


860315 


19.71 


267.00 


77.00 


108.50 


30.10 


30.83 


27.99 




AG 


960513 


19.68 


269.05 


74.12 


105.11 


29.79 


31.09 


30.94 


Doublet #4 


HP 


720525 


19.65 


239.00 


79.00 


113.30 


34.60 


35.05 


23.27 




YK 


911201 


19.62 


235.40 


79.80 


114.60 


34.60 


34.88 


22.22 


Doublet #5 


VR 


610319 


19.73 


154.10 


66.70 


143.00 


44.30 


44.59 


0.35 




HP 


850313 


19.62 


157.00 


65.00 


143.60 


46.30 


46.63 


0.24 


Doublet #6 


HP 


661008 


19.67 


164.00 


50.00 


159.00 


58.80 


61.08 


-5.53 




YK 


750317 


19.67 


163.70 


52.90 


154.90 


56.80 


58.45 


-4.16 


Doublet #7 


HP 


740228 


19.86 


264.00 


58.00 


86.36 


32.52 


41.02 


45.22 




AG 


980330 


19.84 


259.16 


56.32 


84.39 


35.17 


45.44 


45.35 


Doublet #8 


HP 


760206 


19.62 


165.00 


64.00 


140.98 


49.43 


49.49 


2.41 




HP 


850313 


19.62 


157.00 


65.00 


143.60 


46.30 


46.63 


0.24 



Table 6: Event lists of members of clusters within 4° from 4 surface experiments A- 
GASA(AG), Haverah Park(HP), Volcano Ranch(VR) and Yaktuk(YK). 
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Figure 1: Declination distribution of showers of energies above lOEeV and zenith angles 
smaller than 45° of AG, YK and HP. 
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Galactic Latitude ( deg.) Super Galactic Latitude ( deg.) 



Figure 2: Galactic (left) and supergalactic latitude (right) distribution of arrival directions 
of 92 cosmic rays from four ground array experiments. Energy thresholds are 40EeV. Solid 
lines are the sum of expected distributions of each experiment for a uniform distribution. 
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Figure 3: Sky map of 92 events above J^OEeV in equatorial coordinates. Squares - AGASA, 
Triangles - Haverah Park, Circles - Yakutsk and Stars - Volcano Ranch. The region of 
sky observable in each experiment is shown by dashed lines from the top YK, HP, AG and 
VR. 




Figure 4: Sky map of 92 events above J^OEeV in galactic coordinates. The symbols are the 
same as used in figure 3. 
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Figure 5: Frequency distribution of the number of doublets in 1,000,000 simulated data 
sets within several space angles. The arrow mark shows the observed number of doublets. 
Chance probability of the observed number of doublets is calculated by summing the number 
of data sets in the hatched region. 




Figure 6: Arrival directions with error circles of EHECR in supergalactic coordinates. The 
total number of events is 92. 
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Figure 7: Probability of observing doublets and triplets within 1(F from the SG plane with 
1,000,000 simulated data sets for 92 events. 




Figure 8: The location of galaxies within 100 Mpc in G coordinates from the CfA 1995 
catalogue and the clusters within 4°. The squares show triplets and the circles show dou- 
blets. 
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